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We discuss the baryon-number susceptibility \B and related topics which include 
the density fluctuations around the critical point of the chiral transition at fi- 
nite temperature T and baryon density ps ■ Phenomenological implications of the 
density fluctuations near the first-order chiral transition at finite density are also 
discussed. 



1 Introduction 

When exploring a phase transition in any physical system, the study of fluctu- 
ations of physical quantities, especially ones related to the order parameter is 
as important as that of the phase diagram for the system in equilibrium. The 
fluctuations of observablcs arc related with dynamical phenomena including 
the transport properties of the system. The chirally restored and deconfincd 
phase is expected to be created dynamically in the intermediate stage of the 
ultra-relativistic heavy ion collisions and in the early universe. Therefore the 
study of the fluctuations has a great relevance to phenomenology. nn 
In the present report, we discuss the baryon-number susceptibility xb ™ 
and related topics which include the density fluctuations around the critical 
point of the chiral transition at finite temperature T and baryon density pb- 

2 QCD phase diagram in (T, /9s)-plane and the vector coupling 

The lattice simulations of qcdI! suggest that the order and even the existence 
of the phase transition(s) at finite temperature T are largely dependent on the 
number of the active flavors when the physical current quark masses are used: 
For Tji u ~ TTid ~ lOMeV < lOOMeV ~ m s , the phase transition may be weak 
1st order or 2nd order or not exist. The lattice QCD is, unfortunately, still 
not matured enough to predict a definite thing about the phase transition at 
finite baryon density pb(ot chemical potential ijl). 

Low-energy effective modellm and the chiral random-matrix theory] have 
given suggestive pictures of the phase diagram of QCD in the T-fx (or T-ps) 
plane. For example, the NJL modeEl well describes the gross features of the T 
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dependence of the quark condensates of the lightest three quarks as given by 
the lattice QCD, and predict that the chiral transition at fi ^ is of rather 
strong first order at low temperatures T, say, lower than 50 MeV, provided that 
the vector coupling between the quarks as given by g v /2 ■ {qj^q) 2 is absenoij. 
[] As a matter of fact, the strength and even the existence of the 1st order 




sition are strongly dependent on the strength of the vector coupling g v 
the vector term prevents a high-density state. 



The reason why the vector coupling weakens the phase transition and 
postpones the chiral restoration is understood as follows. Thermodynamics 
tells us that when two phases I and II are coexistent, their temperatures Tjjj, 
pressures Fiji and the chemical potentials pi y n are the same; 

T I = r„, Pi = P„, Mi = Mil- (2.1) 



If the phase I (II) is the chirally broken (chirally restored) phase, the last 
equality further tells us something because /iiji at vanishing temperature are 

given by /i; = M 2 + p 2 F , , (i = I, II), where pp i is the Fermi momentum of 

the i-th phase, and M and m are the constituent quark mass and the current 
quark mass that vanishes in the chiral limit. One readily sees that pf { < PF rl , 
accordingly, pb, < Pb u , i-c, the chirally restored phase is in higher density 
than the coexistent broken phase. The vector coupling above gives rise to 
a repulsion proportional to the density squared i.e., gvp 2 B j2 which is bigger 
in the restored phase than in the broken phase. Thus the vector coupling 
weakens and/or postpone the phase transition of the chiral restoration at low 
temperatures. 

Since the vector coupling contribute to the energy repulsively, it also sup- 
press the density fluctuations, or the baryon-number susceptibility, which is 
the main subject of the present report. 



3 Baryon-number susceptibility 

The baryon-number susceptibility \b is the measure of the response of the 
baryon number density on = X)i=i~tv Pi *° infinitesimal changes in the quark 
chemical potentials pn Era: 

N f N f 

XB (T,n) = Ea-](X>) = ((N 2 B ))/VT, (3.1) 

i— 1 i— 1 

c The first order chiral transition in fhhe density direction at low temperatures is also obtained 
in the chiral random-matrix theoryL] 
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where Nb is the baryon- number operator given by Nb = 53i=i N%, with 

Pl = TriV, exp[-0(H - £ f^N^/V = ((N))/V (3.2) 

i— u,d 

the z-th quark- number density, V the volume of the system and (3 — 1/T. 
In the following, we shall confine ourselves to the SU ^(iVj)-symmetric case; 

Mu = Md = Ms = ■•■ = M- 

The baryon- number susceptibility sX Pb ^ is related with the (iso- 
thermal) compressibility of the system □ k t = —Ng 1 (dV/d^)T,N B = 
which tells that if xb is large and so is the density fluctuation, the system 
is easy to compress. 

Another physical meaning of xb is that it is the density-density correlation 
which is nothing but the 0-0 component of the vector-vector correlations or 
fluctuation!]; 

Xb(T, p q ) = (3 J dxS o(0, x), (3.3) 

where S^ v (t,x) = {{j^(t, x)j u (Q, 0)}), with j^x) = q(t,x)'y IM q(t,x) being the 
current operator. Using the fluctuation-dissipation theorem, one has 

X B (T, M9 ) = -limL(0,k), (3.4) 

k — >0 

where L(w,k) is the longitudinal component of the retarded Green's function 
or the response function in the vector channel; 
R tl „(w,k) = F.T.(-i6(t)((\j tt (t,x), jJf,0)]-)))- 

The lattice simulations of QCDO show that xb at /x = is suppressed 
in the low temperature phase, increases with T sharply around the critical 
point of the chiral transition and takes almost the free quark gas value then 
saturates. This behavior may be understood intuitively and roughly as follows: 
In the confined phase at low T, the density fluctuation picks up the Boltzmann 
factor e -Mjv / T with being the nucleon mass, which is much smaller than 
the factor e~ Mq ^ T with M q being the current quark mass (constituent quark 
mass), which factor will be picked up in the deconfined and chirally restored 
(chirally broken) phase. 

Our point here is however that the nature of the chiral transition and 
also presence (or absence) of the vector coupling affect the baryon-number 
susceptibility, especially when ps ^ 0, for which the lattice data is not available 
so fail 
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3.1 Free quark gas 

It is instructive to examine \b (T, fi) in the simple free-quark gas model: 

PB = 2N f N c J Jte(n(T,ri-n(T,ri) (3.5) 

where n(T,n) = l/[exp0(E p - n) + 1] and n(T,fi) = l/[exp (3(E p + fi) + 1] 
with E p = \J M 2 + p 2 , and N c = 3 is the number of the colors. Then one 
readily obtains 



XB(T,fi) = 2N f N c (3 J * |n(l-n)+n(l-ft)} = xi 0) (?», (3.6) 
which is reduced to 



v (Q)r T ^ - f Tl - at a f d P exp{E p /T) 

Xb (T,0) = x B (T) - 4N f N c pJ (3.7) 

at fx = 0. 

If M(T) is decreased as in the chiral restoration, %s increases and reaches 
NfT 2 at M(T) = 0: The enhancement is, however, found to be modest and 
not so large as obtained in the lattice simulations. 



3.2 Model calculation 

To demonstrate the relevance of the nature of the chiral transition and the 
presence of the vector coupling to xb, we perform a calculation with an ef- 
fective model, which is given by adding the vector-coupling terms E3 to the 
Nambu— Jona-Lasinio modeH: 

Nj-l 

£ = q(ij-d-m)q+ y[(9 A a<?) 2 + (qiXalsq) 2 ] 

Nj-l 

"|E mal.q? + (qXal^qf] (3.8) 

a=0 

The realistic value of g v used in the literature i is roughly in the range of 
g v A 2 = 5-9. 

In the self-consistent mean field approximation, the constituent mass (dy- 
namical mass) M,(T, fx q ), the quark condensates ((q~iqi)) and the quark density 
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Pi are all coupled with the vector coupling g v and determined by the following 
equations; 



M i = m i -2g s ((q i q i )), (3.9) 
= -2N C J * {l-n^T.MO-aiCT.fc)}, (3.10) 

Pr = 2N c j -^(ni(T,ii)-ni(T,ii)), (3.11) 

where it is to be noted that the shifted chemical potential p = p — 2g v pi enters 
the distribution functions instead of the naive one, p. 

Simply differentiating these equations with respect to p, one obtains xb(T, p). 
It is noteworthy that when p ^ there arises a coupling between xb and the 
scalar-density susceptibility Xs owing to the non-vanishing "vector-scalar sus- 
ceptibility" Xvs ■ They are defined by 

Xs = =fij d*««(0, x)g(0, x)«(0, 0)q(0, 0))), (3.12) 

Xvs = = / «&((9(0, x) 7 o<z(0, x)g(0, 0)g(0, 0))), (3.13) 

respectively. Xs represents the fluctuation of the order parameter of the chiral 
transition, and is related with the sigma meson propagator. The differentiation 
leads to the coupled equation 



i + 2 9vX f -x<°> )( XB W xf 
4 5s ff.xl ] 1 - 2gsXs 0) ) \ -29sX vs ) \ 2g sX n 



(3.14) 



where x^(T,// 9 ) is the zcro-th order baryon- number susceptibility given be- 
fore, 

N f 

x $i = - m ^i J (2^3 fa 1 - ^ - ^ - ^ 



the zero-th order vector-scalar one and 

N, 



v(0) = 27V V [ -*_[^-(l 



-m-m 



M 2 

+/?-=Hn i (l-n i )+n i (l-n i )}] (3.16) 



the scalar-density susceptibility in the zero-th order. One should note that 
when [i = 0, then Xvs van ishes. O ne can readily obtain Xq ancl Xs hi terms of 
X^°\xvs arL( ^ Xs 1 so we do not write them down to save the space. 



3.3 (A) Finite density case gv = 

To see how the nature of the chiral transition can affect the behavior of xb, 
let us first take a simple case where gv is negligible. 

The essential point lies in the fact that the distribution function n(T, fi) = 
[ e 0( E p-v) _|_ I]- 1 depends on not only explicitly but also implicitly through 
the dynamical mass M(T, /i) in E p = \J M" 1 +p 2 ] hence 

dn , . M dM , , . . 

r-=nl-n---nl-n 3.17 
op E p dfi 

Thus 

JV, 



{n;(l - Uj) - fii{\ - fii)} . (3.18) 



The notable point is the presence of the derivative ^p- . The constituent mass 
M varies with the quark condensate ((qq)) by Eq. (3.9), the order parameter 
of the chiral transition. We saw in §2 that the chiral transition at low temper- 
atures is likely to be of first order in the chemical potential direction. It means 
that the derivative diverges at the critical point at low temperatures, hence 
so does the susceptibility xb- 

3.4 (B) Zero-density case with gv =/= 
Putting fi = into the expressions one getfl 

XB = Xb \> - , (3.19) 

where x^b (T) ^ s ^ ne susceptibility for the free-quark gas. The denominator of 
Xb is essentially the inverse of the propagator of the vector meson in the ring 
approximation at the vanishing four momenta. The above expression shows 
that xb is suppressed by the vector coupling ( g v > 0). This is reasonable 
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at least for a system with a finite //; because the system becomes hard to 
compress when the vector coupling is present, the number fluctuations will be 
suppressed^ Recall also that \b is proportional to the compressibility kt- 

The comparison with the lattice datil shows that the vector coupling is 
rather small in the high temperature phase, ft is interesting that this sup- 
pression of the vector coupling at high temperatures is consistent with the 
observation that-.the screening masses of the vector modes obtained in the lat- 
tice simulations □ almost coincides with 2nT, the lowest screening mass of the 
q-q system in the chiral limit; a similar result is also obtained in the instanton 
approacHa. 

4 Estimate of g v and g s in the lattice QCD 

Boyd et al@ once extracted the effective coupling constants g v and g s from 
the lattice data using the expressions given in the NJL model as given above; 

Xb = X ( b/( 1 +9vX ( b) and Xs = Xs°V( 1 -5 s Xs 0) )- The y concluded that 
when T ~ Tc, g s is much bigger than g v ; g s ~ 4g v . This result is consistent 
with our analysis and the behavior of the screening masses. 

5 Implications to phenomenology 

Kumagai, Miyamura and Sug itat<0 discussed the implications of the baryon 
number susceptibility and the strangeness susceptibility to the observables in 
the relativistic heavy-ion collisions. They argued that in the stopping region 
where the chemical potential is large, large fluctuations of the baryon and the 
strangeness numbers may be a signature of the chiral transition. 

Here we wish to also indicate that the large number fluctuations cause 
those in the scalar channel (the sigma meson channel) at finite density. 

The observability of the possible large densitsz-fluctuations caused as a 
critical phenomenon is discussed by several authored. 

6 Summary and concluding remarks 

We have examined the baryon-number susceptibility \b as an observable which 
reflects the confinement-deconfinement and the chiral phase transitions in hot 
and/or dense hadronic matter. 

The suppression of \b at low temperatures and steep rise around the 
critical temperature as shown in the lattice QCD may be roughly attributed 

d The interactions due to vector fields like u> or neutral p mesons increases the repulsion 
between the constituents of the system. 
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to the confinement-deconfinement transition. Nevertheless, we have shown 
that such a behavior of xb is also affected by the chiral transition. 

Since xb is a measure of the rate of the density fluctuation in the sys- 
tem, the chiral transition at finite chemical potential especially leads to an 
interesting phenomenological consequence to Xb- When the vector coupling is 
small, the chiral transition at low temperatures is of first order in the density 
direction, which implies a divergent behavior of xb, accordingly a huge den- 
sity fluctuations. We have emphasized that such a large enhancement of the 
fluctuation can be also expected for the scalar density fluctuations due to the 
scalar-vector mixing at finite density. Such a large enhancement may leads to 
an enhancement of the sigma-meson production 0. The above phenomena all 
have relevance to experiments to be done in RHIC and LHC. 

We have indicated that the nature of the chiral transition as to the first 
order or not etc is sensitively dependent on the strength of the vector coupling. 
An analysis of the lattice data suggests that the vector coupling is small in 
comparison with the scalar coupling at high temperature. 

The susceptibility xb is nothing but the generalized susceptibility x(w, k) 
at u) = k = 0. One should examine x( w ; ^) m the whole region of lo and k to get 
more information about the vector correlations and the density fluctuations. 

In conclusion, we would like to thank the organizers of this workshop, 
especially, Prof. H. Suganuma, for inviting me to talk on the baryon-numbcr 
susceptibility. 
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